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Using a new quark-field construction algorithm and a large variational basis of operators, we 
extract a highly excited isovector meson spectrum on dynamical anisotropic lattices. We show how 
carefully constructed operators can be used to reliably identify the continuum spin of extracted 
states, overcoming the reduced cubic symmetry of the lattice. Using this method we extract, with 
confidence, excited states, states with exotic quantum numbers (0 + ~, 1 h and 2 + ~) and states of 
high spin, including for the first time in lattice QCD, spin-four states. 



Introduction: The spectroscopy of excited meson 
states is enjoying a renaissance through the observations 
of multiple new states in the charmonium sector. This 
will continue through the forthcoming experimental ef- 
forts at GlueX, BES III and PANDA that will probe the 
spectroscopy of mesons in both the light and charm sec- 
tors. New states demand explanation within QCD and 
may offer insight into the appropriate degrees-of-freedom 
of low energy QCD. A particular example is mesons of 
exotic J PC , those states whose quantum numbers can- 
not be constructed from a quark-antiquark bound state, 
and whose existence may signal the influence of explicit 
gluonic degrees of freedom. 

Lattice QCD provides an ab initio method for the de- 
termination of the hadron spectrum. This approach to 
spectroscopy necessitates methods for measuring the two- 
point correlation functions of field operators with the se- 
lected quantum numbers under investigation. However, 
it has proven difficult to extract precise information from 
lattice QCD about excited states, states of high spin and 
states with exotic J PC . In this letter we will present 
results using a large basis of composite QCD operators 
and a variational analysis method which show that such 
extractions are now possible. 

Access to states of spin-two or higher requires opera- 
tors with spatially separated quark fields. To facilitate 
this kind of construction, a new quark-field construc- 
tion algorithm, called "distillation" , was developed [l| 
recently which enables efficient calculations of a broad 
range of hadron correlation functions, including those 
with spatially separated quark fields. 

In Euclidean space, excited-state contributions to cor- 
relation functions decay faster than the ground-state and 
at large times are swamped by the larger signals of lower 
states. In improving our ability to extract excited states, 
better temporal resolution of correlation functions proves 
extremely helpful. An anisotropic lattice, where the tem- 
poral direction is discretized with a finer grid spacing 
than its spatial counterparts, is one means to provide this 
resolution while avoiding the increase in computational 
cost that would come from reducing the spacing in all 



directions. To this end, a large-scale program has been 
initiated to generate dynamical anisotropic gauge fields 
with two light clover quarks and one strange quark 0, 0] . 

In this work, these anisotropic lattices are combined 
with the distillation technique for the construction of 
quark-antiquark operators with multiple derivative inser- 
tions. Only the connected Wick contractions arc com- 
puted, giving access to isovector quantum numbers. For 
this first investigation, the three-flavor degenerate-quark- 
mass dataset is used (m^ = rriK = m v m 800 MeV) , with 
lattice spacings a s ~ 0.12 fm, a^ 1 ~ 5.6 GeV and a spa- 
tial lattice extent of ~ 2fm. We will argue later that 
using a relatively large quark mass in this first study re- 
duces complications due to mesons being able to decay 
into multi-meson states. 

Spin on a cubic lattice: Lattice QCD computations 
consider the theory discrctiscd on a four-dimensional Eu- 
clidean hypercubic grid. The reduced three-dimensional 
rotational symmetry with respect to the continuum intro- 
duces complications when one wishes to study particles 
of a particular spin, since spin no longer identifies irre- 
ducible representations of the cubic symmetry group 
There are five single-cover lattice irreps for each parity 
and charge-conjugation: Ai,Ti,T2,E,A%. The distribu- 
tion of the various M components of a spin- J meson into 
the lattice irreps is known as subduction, the result of 
which is displayed in TablcQ] In the continuum limit, the 
full rotational symmetry is restored and the components 
subduced into different irreps will be degenerate, whereas 
at finite lattice spacing they will be split by an amount 
scaling with at least one power of the lattice spacing, a s . 

This suggests a simple method to assign continuum 
spins by attempting to identify degeneracies across lattice 
irreps compatible with the subduction patterns in Table 
HI Unfortunately the empirical meson spectrum shows a 
number of approximate degeneracies that may be con- 
fused with those originating through discretisation. As 
an example consider the Xc0,i,2 states in charmonium, 
split only by a small spin-orbit force. These states would 
appear in a lattice computation as a single state in each 
of A^ + , T x ++ , T% + and E ++ and could easily be mistaken 
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J irreps 

Ai(l) 

1 Ti(3) 

2 T 2 (3)ffi£(2) 

3 ri(3)©T 2 (3)®A 2 (l) 

4 J 4i(l)eT 1 (3)ffiT 2 (3)®£(2) 

TABLE I: Continuum spins subduced into lattice irreps 
A(dim). 

with a single J = 4 state split by discretisation effects. 
In the high lying part of the calculated spectrum, shown 
in Figures [T] and [2] we observe considerable degeneracy 
that renders spin-identification by this method virtually 
impossible. In this letter we consider using the addi- 
tional information embedded in the overlaps of states on 
to carefully constructed operators at zero momentum. 

Meson operators: By using a circular basis for both 
spatial derivatives and the three-vector-like gamma ma- 
trices (7i), we can utilise the SO (3) Clebsch-Gordan co- 
efficients to construct continuum operators of definite 
spin. For example, with one derivative and a vector 
gamma matrix we can construct operators of overall spin 
J = 0,1,2: 

(r x D™ =1 ) J - M = (l,mi; l,m 2 | J,M) $T mi *D m J 

where repeated m indices are summed. In the distillation 
framework, the fermion fields, ip are smeared using a low- 
rank filtering operator. 

In the case of two derivatives we couple into a definite 
spin before coupling to the gamma matrix: 

(T x D [ ? d ) j > m = (l,m 3 ; J D ,m D \j,M) 

x (l,mi; l,m 2 | Jd,tod) x ^r„ i3 D mi D m2 tp. 

For three derivatives combining the outermost derivatives 
together first ensures definite charge-conjugation: 

(T x D l $ StJD ) J ' M = (1, m 4 ; J D ,m D \j, M) 

x (l,m 2 ;Ji3,m 13 \j D ,m D ) x (l,mi\ 1, m 3 | Ji3,mi 3 ) 

This scheme can be extended to any desired number 
of covariant derivatives, which in practical computa- 
tions are replaced by gauge-covariant finite differences. 
The gauge links appearing in these differences are stout- 
smeared to reduce UV fluctuations. To be of any real 
use in lattice calculations these operators of definite con- 
tinuum spin, J, must be subduced into the irreducible 
representations of the cubic lattice rotation group (A = 
{Ai, Ti, T 2 , E, A 2 }). Noting that each class of operator is 
closed under rotations, the subductions can be performed 
using known linear combinations of the M components 



for each J: 

52s$ (r x d^Y m = J2 s A,f° J,M 

M M 

where A is the "row" of the irrep. Note that, although 
O^j can have an overlap with all spins contained within 
A, it still carries the memory of the J from which it was 
subduced, a feature we exploit below. 

Spectral analysis: For each lattice irrep A PC the full 
matrix of correlators Cij(t) was computed with equiva- 
lent rows (A) averaged over. The dimension of the matrix 
is therefore equal to the number of operators constructed 
in that irrep. 

The correlation matrix can be described by a spec- 

tral decomposition Cij(t) = Xm 2m 3 e~ mnt (we only 
consider zero momentum), where Zf = (O|0j|n) en- 
codes the overlap of state n on to operator Oi. An 
optimal method (in the variational sense @, @]) to ex- 
tract mass and Z information from the matrix of corre- 
lators is by solution of a generalised eigenvalue problem, 
Cij(t)vj = X n (t 7 t )Cij(to)Vj where the eigenvectors v n 
are related to the Z by Zf = ^/2rn^e m ^°/ 2 vf*C i:j (t^i. 
Our implementation of this approach is described in [7[. 

The extracted spectrum across lattice irreps, including 
all operators with up to three derivatives, is shown in 
Figures Q] and [21 We give the number of operators in 
each irrep and the color coding indicates continuum spin 
assignment suggested by a method we now describe. 

Our particular choice of operator construction of- 
fers us a method to identify the continuum spin of a 
state. We take advantage of the fact that, at the lat- 
tice spacing we work, we expect lattice operators act- 
ing on extended objects such as mesons to behave in 
a manner reasonably close to the full rotational sym- 
metry. In the continuum our operators are of definite 
spin such that (0\O J < M \ J', M') = Z^8j t j.8 M ,M' and so 
(0\O [ £ x \J',M) = SjgzWSjj,. Z^ is a single number 
of dynamical origin describing the overlap of the state of 
spin J on to the operator used. We form a correlator in 
a given irrep A and average over equivalent rows, A, 

dim(A) X! ^AA = dim(A) ^ (° I ^A.A^A.A I °) ' 
A A 

Inserting a complete set of meson states between the op- 
erators and using the fact that the subduction coeffi- 
cients form an orthonormal matrix, J^m <^a X ^A'^X* = 
<5a,a"5a,A', we obtain terms proportional to Z^*Z^; 
these terms do not depend upon which A we have sub- 
duced into. Hence for example a J = 3 meson created 
by a [J = 3] operator will have the same Z value in 
each of the Ai,T\,T<i irreps. Since this derivation uses 
smoothed, semi-classical fields it is valid in the contin- 
uum limit and at finite lattice spacing we expect there 
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FIG. 1: Extracted spectrum of states in the PC = + — ,++ 
channels displayed by lattice irrep. The number of operators 
in each irrep is given below the irrep label. All masses scaled 
by the Q baryon mass as extracted on this lattice 0]. Boxes 
represent the extracted mass and one sigma statistical uncer- 
tainties. Color coding indicates continuum spin identification. 
Orange boxes have well determined masses but undetermined 
spin. Grey boxes have masses that are not well determined 
by the variational fitting method. States with exotic quantum 
numbers + ~ and 2 + ~ are highlighted. 
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FIG. 2: As previous but for PC = — h, 
exotic 1 _+ is highlighted. 



The lowest lying 



to be small deviations from equality due to discretisation 
effects. 

We take advantage of these properties to identify the 
spin of the extracted states in the following way. Firstly, 
we consider the relative magnitudes of the extracted Z 
values for various states. Figure [3] shows that for the 
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FIG. 3: Overlaps, of a selection of operators on to states in 
each lattice irrep, A~~. Z's are normalised so that the largest 
value across all states is equal to 1. Lighter area at the head 
of each bar represents the one sigma statistical uncertainly. 



J mesons, each state has large overlap only onto op- 
erators of a single spin. The second stage of the iden- 
tification requires us to match states in different irrcps 
and compare their Z values with common operators sub- 
duced across irreps. As shown in Figure 01 these values 
agree well. Any small discrepancy could be attributed 
to two causes: discretisation errors from the use of sim- 
ple central-difference operators to represent derivatives 
or the effect of renormalisation. These operators act on 
smoothed gluonic and quarks fields and this eliminates 
fluctuations at the cut-off scale so the latter effects will 
most likely be very small. 

Using this method we have extracted a large number of 
states with all possible PC combinations and confidently 
identified the spin of these states; the spectrum is shown 
in Figures [1] and [2j As well as extracting many excited 
states, we have for the first time identified states with 
spin four: 4 ++ , 4 h and 4 . We have extracted states 
with exotic quantum numbers (0 -1 , 1 ^ and 2^ ) and 
these are highlighted in the figures. The presence of these 
exotics likely points to the influence of explicit gluonic 
degrees of freedom. 

Figure [3] shows that the third excited vector state 
(m/mo ~ 1.35, marked with f) has a qualitatively differ- 
ent pattern of Z values compared to the lighter vectors. 
Notably, there is large overlap with the (tt x D^fL^^ 1 
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FIG. 4: A selection of Z values across irreps A for states suspected of being J = 2, 3, 4 respectively. 



operator. -Dji, corresponds to the commutator of two 
covariant derivatives which vanishes in the absence of a 
gluonic field. This commutator is proportional to the 
field strength tensor and so the significant overlap hints 
at a gluonic component. This suggests an identification 
of this state as a crypto-exotic vector hybrid, although 
the non-zero overlap onto ip^fiip suggests some mixing 
with a conventional vector state. 

Two-meson states: We might expect to observe an 
abundance of two-meson states above 2/71^ ~ 0.85mn, 
but such states are not apparent in our extracted spec- 
trum. This is most clearly seen in the A" channel 
where the lightest state extracted is a J = 4 state above 
1.5mn, while a pseudoscalar- vector state with the mini- 
mum relative momentum allowed in our finite box would 
be expected close to 1.2ma- The operators used in this 
study featured only a single ip, ip field pair and so do 
not have overlap onto quark Fock states higher than qq. 
QCD dynamics can act to mix qq Fock states with two- 
meson basis states to form mesonic eigenstates. This 
mixing is expected to be significant when a discrete lat- 
tice two-meson state is degenerate with a "single meson" 
to within that meson's continuum decay width. At this 
relatively heavy quark mass, we expect low-lying reso- 
nances to have small widths due to reduced phase-space 
for their decay and hence for there to be only small mix- 
ing with two-meson states, perhaps explaining our lack 
of observation of such states. A calculation similar to 
the one reported herein has been carried out on a lat- 
tice of spatial extent ~ 2.4 fm. The extracted spectrum 
is found to be identical within statistical fluctuations to 
that presented here. This is more evidence that we are 
not seeing two-meson states since their allowed relative 
momentum, and hence their energy levels, would have 
changed significantly. These issues can be properly in- 
vestigated by including in the variational basis operators 
featuring a product of two fermion bilinears, expected to 
have good overlap on to two-meson states. This work is 
underway 



Summary: We have demonstrated a lattice QCD op- 
erator construction that enables the identification of con- 
tinuum spin with some confidence. Using distillation 
technology to construct the correlators, and a variational 
analysis to study them, we have extracted an excited 
state spectrum featuring well-determined states with ex- 
otic quantum numbers and, for the first time, states of 
spin-4. 

It is notable that our extracted spectrum has both fea- 
tures of the n 2S+1 Lj state assignment of bound-state 
quark models and also states that do not seem to lie 
within that classification. We believe that this study is 
seeing a full spectrum of QCD mesons which includes 
exotic and non-exotic hybrid mesons Q. 
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